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We study scaling of the ground-state fidelity in neighborhoods of quantum critical 

points in a model of interacting spinful fermions - a BCS-like model. Due to the exact 
diagonalizability of the model, in one and higher dimensions, scaling of the ground- 
state fidelity can be analyzed numerically with great accuracy, not only for small 
systems but also for macroscopic ones, together with the crossover region between 
them. Additionally, in one-dimensional case we have been able to derive a number 
of analytical formulae for fidelity and show that they accurately fit our numerical 
results; these results are reported in the article. Besides regular critical points and 
their neighborhoods, where well-known scaling laws are obeyed, there is the multi- 
critical point and critical points in its proximity where anomalous scaling behavior is 
found. We consider also scaling of fidelity in neighborhoods of critical points where 
CN , fidelity oscillates strongly as the system size or the chemical potential is varied. Our 

jfT results for a one-dimensional version of a BCS-like model are compared with those 

l/-) ■ obtained by Rams and Damski [IT] in similar studies of a quantum spin chain - an 

I/"") , anisotropic XY model in transverse magnetic field. 

^ : 1 Introduction 

In the last decade quantum phase transitions and quantum critical phenomena continue to 
be a subject of great interest, vigorously studied in condensed matter physics. Both, ex- 
perimental and theoretical developments point out to the crucial role that quantum phase 
transitions play in physics of high-T c superconductors, rare-earth magnetic systems, heavy- 
fermion systems or two-dimensional electrons liquids exibiting fractional quantum Hall effect 
[I], [2]. The so called classical, thermal phase transitions originate from thermal fluctuations 
and are mathematically manifested as singularities in temperature and other thermodynamic 
parameters of various thermodynamic functions, and such correlation quantities like the cor- 
relation length, at nonzero temperatures. In contrast, quantum phase transitions originate 
from purely quantum fluctuations and are mathematically manifested as singularities in sys- 
tem parameters of the ground-state energy density, which is also the zero-temperature limit 
of the internal energy density. Naturally, singularities of thermodynamic functions appear 
only in the thermodynamic limit. The importance of quantum phase transitions for physics 
and the related wide interest in such transitions stems from the fact that, while a quantum 
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phase transition is exhibited by ground states, hence often termed a zero-temperature phe- 
nomenon, its existence in a system exerts a great impact on the behavior of that system also 
at nonzero temperatures. 

Typically, theoretical studies of quantum phase transitions can be conducted much along 
the same lines as in the case of thermal phase transitions: one considers possibly simple mod- 
els, studies the eigenvalue problem of Hamiltonians, excitation gaps, constructs local order 
parameters, determines broken symmetries, calculates the correlation functions. However 
in distinction to thermal phase transitions, completely new approaches to the subject have 
been put forward, which are based on ideas of the quantum-information science. The cen- 
tral object of quantum-information approaches to many-body systems is the ground state; 
there is no need to consider order parameters and broken symmetries, which makes these 
approaches particularly useful for quantum phase transitions, which are not characterized in 
such terms, like topological phase transitions. One of the quantum-information approaches, 
the entanglement approach, with longer than a decade history, exploits the notion of en- 
tanglement of quantum states in the context of its applications to many-body systems, for 
review see [3J; this approach is out of the scope of our article. Another approach, the fidelity 
approach, put forward by Zanardi and Paunkovic [I], is a relatively new idea of using the 
quantum fidelity of ground states as a probe of quantum phase transitions. A comprehensive 
exposition of the fidelity approach can be found in the review by Gu [5]. 

In what follows we shall study continuous quantum phase transitions, and it is the fidelity 
approach that will be used for that purpose. Let A be a vector whose components are those 
parameters of the considered system's Hamiltonian that drive a quantum phase transition, 
and e - a unit vector in the space of those parameters. Then, on varying parameter S the 
vectors A + 5e scan a neighborhood of A along direction e. For given A and 5, the ground- 
state fidelity at A in direction e, ^(A, 5), is the absolute value of the overlap of the ground 
states | A ± 5e) at the points A ± Se, 

3 r e (\,S) = \(\-Se\\ + Se)\. (1) 

A list of general, system independent, properties of ^(A, 8) can be found in [5]. The 
transition point of a continuous quantum phase transition, i.e. the quantum critical point, 
denoted A c , is characterized by the power-law divergence of the correlation length £(A), 
as the quantum critical point is approached: £(A) ~ |A — \ c \~ u ■ Alternatively, A c can 
be defined as the point where the gap between the ground-state energy and the energy of 
the lowest excited state vanishes. In reference [I] it was demonstrated that the quantum 
critical point A c can be identified as the minimum of fidelity, as A is varied. However, the 
fidelity approach seeks an answer to a more general question: does the behavior of quantum 
fidelity in a neighborhood of a quantum critical point encode not only the location of that 
point but also some universal information about the underlying quantum phase transition? 
First results, pointing towards a positive answer to the raised question by providing some 
finite-size critical scaling of geometric tensors, have been obtained by Venuti and Zanardi 

According to finite-size scaling theories, the properties of a system are close to those at 
the thermodynamic limit, we say the system is macroscopic, if the linear size of the system, 
L, is much greater than £(A). In the opposite limit we observe the so called "small-system" 
properties. Concerning finite-size scaling properties of fidelity 3^ (A, S), it is expected that 
the characteristic length of the system, that differentiates between small and macroscopic 
system, is given by £ e (A, S), which is the smaller of the two correlation lengths £(A ± Se). In 
other words, the crossover between small-system and macroscopic-system properties occurs, 



when the effective linear size of a system, L/£ e (A, £), satisfies the crossover condition: 

L/£ e (A,<S)~l. (2) 

There are numerous papers devoted to critical scaling of small-system fidelity, see [5] , [6] , [7] , 
and references quoted there. Typically, whether A is close to or away from A c , small-system 
fidelity can be Taylor-expanded in 5, 

y.(A,«y) = l-yx«(A) + ..., (3) 

where the first order term vanishes because of the symmetry of fidelity in S at zero. The 
coefficient of the second order term, x e (A), is known as the fidelity susceptibility. One expects 
some universal scaling properties of Xe(A), provided A is sufficiently close to a quantum 
critical point A c , where the correlation length diverges: £(A C ± 8e) ~ \8\~ u - Fairly general, 
model- independent, arguments provide us with finite-size scaling of the fidelity susceptibility 
at A c 0: 

Xe (A c ) ~ L 2 I\ (4) 

that is, in the small-system regime 

-ln%(\ c ,5)~5 2 L 2 /». (5) 

In the macroscopic-system regime, quantum phase transitions have been studied by means of 
the so called fidelity per site, a quantity whose logarithm is equal to iV" 1 In 3^ (A, S), where 
N = L d is the number of sites in a <i-dimensional system [8], [S]. In this regime expansion 
([3]) is not valid near critical points [S]-[TT]. However, critical scaling of a macroscopic-system 
fidelity has been considered only very recently by Rams and Damski [TD],[TTJ. These authors 
found that, while for small-systems the fidelity scaling is totaly insensitive to the way the 
critical point A c is approached by the points A ± 5e (i.e. for instance, whether they are 
located on one side of the critical point or on the opposite sides), in the case of macroscopic- 
system the way of approaching the critical point matters. To make this explicit, Rams and 
Damski substituted A c + c5e for A. By choosing the value of the parameter c, the above 
mentioned location of the two points A c + c5e ± 5e with respect to the critical point can 
be controlled. If \c\ > 1 or \c\ < 1, then both points are located on one side or on opposite 
sides of A c , respectively. If \c\ = 1, then one of the points coincides with A c . Now, the above 
mentioned independence of the small-system-fidelity scaling on the way the critical point A c 
is approached by the points A ± Se can be expressed as follows: 

-ln3' e (\ c + c8e,5)~5 2 L 2/v . (6) 

In contrast to the small-system case, the fidelity scaling law for macroscopic systems, derived 
by Rams and Damski [TU], makes the dependence on parameter c explicit, and reads 

-ln? e (\ c + c5e,5) ~ \5\ du NA e {c), (7) 

where A e (c) is the scaling function. Scaling law (^ is expected to hold provided the ther- 
modynamic limit of N~ 1 ln3 r e (X,5) does exist, there is only one characteristic length scale 
given by £, and du < 2 [TOj.pTj. 

The testing ground for new developments in theory of quantum phase transitions con- 
sists mainly of two paradigmatic quantum spin chains, an Ising chain and an (isotropic or 
anisotropic) XY chain, both in a transverse magnetic field. In this respect, the fidelity ap- 
proach is no exception as can be seen in references [5] - [TT]. This stems from the exact 



diagonalizability of both models, which makes possible some analytical analysis. In those 
models where such an analysis is not feasible, only rather small systems are studied nu- 
merically, with no perspective to reach a macroscopic regime [5]. Due to Jordan- Wigner 
transformation the two mentioned quantum spin chains are equivalent to one- dimensional 
spinless-fermion lattice gases that are exactly diagonalizable for a specific boundary condi- 
tion. Beyond one dimension those models are not exactly diagonalizable, hence less useful 
for testing theories of quantum phase transitions. 

To test the fidelity approach, in particular critical scaling of fidelity, from small to macro- 
scopic systems, we consider spinful fermion models that are exactly diagonalizable in any 
dimension, for periodic, antiperiodic or twisted boundary conditions. These models originate 
from the two-dimensional model of d-wave superconductivity proposed by Sachdev [T2|, [T]. 
Naturally, the restriction to one-dimension is of special merit, since besides accurate numeri- 
cal results, we have been able to derive analytical results for two-point correlations functions, 
the correlation length, and fidelities in small- and macroscopic-system regimes. These results 
constitute the contents of our paper. Results for two-dimensional systems will be presented 
elsewhere. 

The paper is organized as follows. In section [2] we present our model and basic facts rele- 
vant for subsequent discussion. Then, in 6 subsections of section [3], we study critical scaling 
of fidelity, in small- and macroscopic-system regimes, and the crossover region between them, 
in vicinities of different kinds of quantum critical points of the model. The obtained results 
are summarized in section HI 

2 The model, its quantum critical points, ground- 
state two-point correlation function and fidelity 

We consider a d- dimensional spinful fermion model, given by the Hamiltonian, 

H = ~2 Yl ( a la a l+e„a + h.C. - ^Ol,^ - - ^ 4 (4W+e*4 ~ °U°i+e*,t) + hX ' ( 8 ) 

l,i,cr l,i 

where a l(T , ai t0 stand for creation and annihilation operators of an electron with spin pro- 
jection a =t, 4- i n a state localized at site I of a <i-dimensional hypercubic lattice with some 
boundary conditions, respectively, and e iy with % = l,...,d are unit vectors whose m-th 
component is 8^ m . Both summations in (|S]) are over pairs of nearest neighbors, with each 
pair counted once. The real and positive parameter t is the nearest-neighbor hoping in- 
tensity, /i - the chemical potential, and Jj - the coupling constant in direction d of the 
gauge-symmetry breaking interaction, in general complex. Naturally, we can express the 
parameters \x and J; in units of t, while the lengths of the underlying lattice in units of 
the lattice constant, preserving the original notation. We emphasize that in distinction to 
[12], PP, where Hamiltonian flS} was derived, here the parameters \x and Jj are independent. 
Moreover, we impose on the system the antiperiodic boundary conditions, for which the grid 
of wave vectors k (quasimomenta) is given as follows: a component ki of k assumes the 
values ki = (2n + 1 — L)ir/L, with n — 0, 1, . . . , L — 1, where L is the size of the underlying 
lattice in direction e i; and we choose L divisible by 4. The values of ki are distributed sym- 
metrically about zero, the closest to zero points are ±tt/L and the closest to ±7r/2 points are 
±7r/2±7r/L. Then, after passing from the site-localized to the plane- wave basis, Hamiltonian 



assumes the form 



H = J2 £ kcl,aCk,a ~ J2 [ Yl Ji C ° S ki ) C lt C -fe4 + hX -' 
k,a k \ i / 



where 



Ek = 2_, cos h ~ /•*• 



(9) 



(10) 



Formally, Hamiltonian §§§ differs from the well-known BCS Hamiltonian of s-wave super- 
conductivity by the presence of cos ki factors in the gauge-symmetry breaking term. Such 
Hamiltonians can readily be diagonalized by means of the Bogoliubov transformation. The 
dispersion relation of quasi-particles reads 



Ek + 2_^ ( £ k ~ Ek) 



(11) 



where 



Ek 



\ e k 



/ Ji cos ki 



(12) 



The Hamiltonian (jHJ) preserves parity; therefore without any loss of generality we can restrict 
the state-space to the subspace of even number of particles (electrons). In this subspace, the 
state \0) qp of an unspecified (but even) number of electrons, defined by 



I°)<zp= ]1k + v1, 



t u -k,i 



)|0), 



(13) 



where |0) is the electron vacuum, with Uk real and positive, 



iik 




(14) 



and, in general, complex Vk, 



\ v k\ 




j , arg v k = arg ^ J f cos k t 



(15) 



is the eigenstate of (|HJ) to the lowest eigenenergy, J2k ( £ k ~ Ek), he. the ground state. This 
state is the vacuum of elementary excitations (quasi-particles), whose energies, relative to 
the ground-state energy, are Ek- As long as our system is finite, the energies Ek are the 
excitation gaps, since Ek > for all values of k. However, in the thermodynamic limit, 
L — > oo, the excitation gaps Ek in the spectrum of quasi-particles may close at special 
values of wave vector k, denoted k c , and at special values of the chemical potential /i and 
coupling constants Jj. 

One can define several two-point correlation functions for the considered system, which 
however are simply related with each other. Each one can be used to calculate the correlation 
length. We choose the two-point function given by the matrix elements of the ground-state 
one-body reduced density operator, gp (0|ao i0 .a ri(J |0)gp. Explicitly, 



'//< 



(0|a* a P , ff |0> 



qp 



1 
2L d 



£k 



y — exp ikr. 



(16) 



Multicritical point 




Critical lines 



Figure 1: Critical lines in the (//, J) plane; the versors e show directions in which neighbor- 
hoods of critical points are scanned by the points A c + (c ± l)Se. 



Using the invariance of functions £& and Ek with respect to reflections in coordinate axes, 
and then taking the thermodynamic limit 



1 r d 

lim qp {0\al a r ^\0} qp = 3(r) = -— / dfc-^-TTcosfyr,-, 



(17) 



where r = \r\. 

Let \0) qp and \0) qp be two ground states, the first one for parameters /i, Jj, and the 
functions Ek, Ek, the second one for parameters p,, Ji, and the functions ik, Ek- As a result 
of the product structure of the ground states, the quantum fidelity for these states has also 
a product structure, 



\<iv 



(0|0) 



<n>\ 



\ | (ukUk + \vk\\vk\ expi(&igv k - arg^A,)) 

k 



Further considerations will be restricted to the one- dimensional case, that is in (^ we 
set the coupling constants Jj = 0, for i > 1, J\ = J -- a real number, with J independent 
of /x, and the wave vector k reduces to the wave number k. Consequently, Hamiltonian (jSj) 
is parameterized by the pairs (/i, J). After using (IH1) . f|T5|) the fidelity of two ground states 
(TT3l) assumes the form 



i'//' 



(0|0) 



(/p i 



I/(*), /(*) 



fc>0 



1 I ^fc^fc + J J cos 2 k 

2 \ -Efc-^fc 



(19) 



where the restriction of the product to k > is a consequence of the symmetry in k at zero 
of the factors in product (ITHj) . 

The critical points of the above specified one- dimensional version of (JSJ) are located at 
the line ji = for A; c = ±7r/2, and at the line J = 0, if — 1 < fj, < 1, for k c — arccos(/i). 
We note on passing that these lines are symmetry lines of Hamiltonian (JSJ): at // = line 
if is hole-particle invariant and at J = line it is gauge invariant. We do not point out 
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order parameters able to distinguish ground-state phases, since in the fidelity approach to 
quantum phase transitions they do not play any role (which is one of the advantages of fidelity 
approach). Among the mentioned critical points one can distinguish the multicritical point, 
/i = and J = 0, and two critical end points, /i = ±1 and J = 0, see FigJTJ 

Since for the considerations in the sequel the knowledge of the dependence of the cor- 
relation length on n and J plays a key role we present a summary of our results obtained 
in the one-dimensional case [13]. In a neighborhood of the line of critical points, /i = 0, 
the large-distance asymptotics of the integral in (TT7|) has been determined in [13] by means 
of coupled analytical and numerical arguments. It is convenient to introduce an auxiliary 
function S(?~) : 

^/ n 1 H ,, cos kr , , 

S(r) = -7T- / dk-—, (20) 



2tt 



El 



in terms of which 



S(r) 



S(r + l) + §(r-l) -/i§(r 



(21) 



The large-distance asymptotics of S(r) reads 



§(r) = e(-l) r cos 



where the distance-independent factor C is 



2tt 3 



1/2 



r-^e-'K, 



(22) 



(l_ /3 2 ) -l/4 [(l + M) 2 + J 2 ] -l/ 2 



(23) 



with 







/' 



1 + J 2 



1 + J 2 



arccos 



and 



^* 



(-aO 

i + j 2 ' 

j 2 



j 2 



1/2 



| J|/i 



(24) 



(25) 



(26) 



In formulae (J23|) -( p6|) . the chemical potential is taken positive, ji > 0; they hold for /i/(l + J 2 ) 
sufficiently small (less than 10~ 2 in our calculations). Apparently, the quantity £ that follows 
from formula (1221 for the subsidiary function S(r) is the correlation length, as defined by 
the large-distance asymptotic behavior of the correlation function S(^), in a neighborhood 
of critical line /i = 0. As fi — > the correlation length diverges with the critical index v — 1. 

Unfortunately, we have not been able to find analytic formulae for the correlation length 
in neighborhoods of critical points located at the line J = in the interval — 1 < fj, < 1. 
Numerical calculations show that for < \i < 1, £ = C(/i)|J| -1 , with some function C(aO) 
that is again v — 1. However, at the end critical points, /i = ±1 and J = 0, numerical 
calculations give £ ~ | J| _1 / 2 , i.e. v = 1/2. 

Let us adapt the general notation introduced in section [Q to the considered here model. 
As the location of critical points is uniquely determined by pairs (/i, J), we set A = (/i, J), 
hence |A) = \0) qp . Then, in formula (TT9|) for fidelity, the functions £&, Ek, given by (flOl) and 
(TT2I) . respectively, are calculated at A c + (c — l)£e, while ejt and Ek - at A c + (c + l)£e. 
Finally, we set | 9P (0|0) 9P | = ^(Ac, 5). 



In what follows, we shall study numerically the sum 

J2(-lnf(k)) = -\n%(\ c ,S), (27) 

fc>0 

as a function of parameter 8 for fixed system size L, or vice versa, in neighborhoods of 
various critical points, in small- and macroscopic-system regimes. In all the considered cases 
the function ln/(/c) is either continuous in the whole interval [0,7r] or it has an integrable 
singularity at some k (a discontinuity or a logarithmic divergence). Therefore, in all the 
considered cases the limit L — > oo of the Riemann sum corresponding to ( 1271) does exist, 

lim -L- 1 In 5 e (A c , 8) = -*- f dk (- In f(k)) . (28) 

L-^oo Z7T Jo 

Consequently, for given sufficiently small \8\ and sufficiently large L 

-\n? e (\ c ,5) « A [ dk (-lnj(fc)) (29) 



2-k J0 

approximately, that is in a macroscopic-system regime — In 3e(A c , 8) scales with the system 
size as L. 

3 Critical scaling of the ground-state fidelity 

Any study of critical behavior involves specifying a critical region, that is a critical point 
and its neighborhood. The quantum critical points considered in our paper are displayed in 
Fig. [TJ As for neighborhoods, we have chosen line neighborhoods, each one specified by a 
unit vector e and a range of parameter 8, which are scanned by vectors A c + (c — l)8e and 
A c + (c+ l)8e on varying 8 (see section [1]). It might be expected, and it is indeed the case as 
our studies show, that there is several kinds of specific critical scaling of fidelity, depending 
on the choice of the critical region. Therefore, our results concerning critical scaling of the 
ground-state fidelity are presented in six subsections, each one labeled by the specific location 
of the considered critical points A c and direction e of the line neighborhood. 

In the sequel, we typically calculate — In ^(Ac, 8) as a function of 8, keeping the linear size 
L fixed, for a number of L; in fact — ln^^Ac, 8) depends on \8\, since 3 e (\ c , 8) = ^(Ae, —8). 
However, sometimes it turned out to be instructive to fix 8 and vary L. According to our 
results concerning the correlation length, summarized in section [21 in particular formula 
(12"oT) . for a given linear size of the system, we enter the regime of small system by decreasing 
sufficiently \8\ or the regime of macroscopic system - by increasing it sufficiently. Then, 
we verify the scaling with L and with 8, and identify small-system, macroscopic-system and 
crossover regimes. In many cases we provide analytic approximate formulae for — In ^(Ae, 8) 
in small- and macroscopic-system regimes, and compare them with numerical data. In 
discussing critical scaling properties of fidelity, it is always instructive to study the behavior 
of the function f(k) of formula f|T9l) in neighborhoods of wave numbers k c for which the 
excitation gap Ek closes. 

Finally, let us note that there is no critical scaling of interest for critical points \y\ < 1 at 
the line J = with neighborhoods along the //-axis (e = (1, 0)). In this case the considered 
system is free, its ground state is a Fermi sea and fidelity 3 r (i )0 )((/i, 0), 8) can assume only 
two values, zero or one. The fidelity vanishes for sufficiently large systems, unless both of 
the points (//, 0) + (c ± 1)8(1, 0) are located on one side of the |/x| < 1 interval. 
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Figure 2: (Color online). The plots of / (left panel) and —In/ (right panel) versus tt/2 — k 
for three values of c: c = - red solid line, c = 1 - green dashed line, c = 2 - magenta dotted 
line. Right panel is in doubly logarithmic scale, and the blue dotted straight line indicates 
the — ln3 r (i j0 )((0, J), 5) ~ k~ 2 behavior, and is given by formula (132]) . 

3.1 Across the fi—0 critical line: A c =(0, J), J away from 0, e=(l,0) 

At critical points A c = (0, J), independently of J, the excitation gap E k closes at wave 
number k c = tt/2. In line neighborhoods in direction e = (1,0) of those critical points the 
function f(k) assumes the form 



/(*) = 2 < 



(cos k — n) (cos k — p) + J 2 cos 2 k 



a/[(cos k — /i) 2 + J 2 cos 2 k] [(cos k — /i) 2 + J 2 cos 2 k] 



(30) 



where we set 



f 



{c-l)S, fi = (c+l)6, 



(31) 



with c > throughout this section. 

In this section we consider J being away from zero, say \J\ > 0.1. In Fig. [2] it is shown 
that f(k) exhibits large variations in a vicinity of k c = tt/2 and in this region it is sensitive to 
the value of c; for c = 1, f(k) exhibits a discontinuity at tt/2. Far away from tt/2 the function 
f(k) is insensitive to c. The plot of — ln/(fc) in doubly logarithmic scale (right panel of Fig. 
[2]) reveals the —In f(k) ~ k~ 2 behavior {k = tt/2 — k), sufficiently far away from tt/2. Since 
for given § and sufficiently small system size L, the grid of wave numbers does not sample 
well the vicinity of tt/2, where f(k) varies strongly and is sensitive to c (the closest to tt/2 
points of the grid, tt/2±tt/L are beyond this vicinity), — ln3 r ( 10 )((0, J), 5)) does not depend 
on c and scales with the system size as L 2 . This is the small-system behavior. In Fig. [31 
where L is fixed, this behavior is well visible for sufficiently small 5. To get an approximate 
formula for — ln3 r ( 10 )((0, J), 5)) in the small-system regime, in formula (130]) we approximate 
cos k by k and expand in 5 at zero, keeping only the lowest order term: 



- In /(tt/2 -k) 



J 5 



J 2 K 



(32) 



and then we sum contributions from the two closest to tt/2 wave numbers k = ±tt/L: 



-ln% )((0,J),5) 



2_ J 2 

^(1 + J 2 ) 



;(8Lf 



(33) 



Formula (I33p approximates well — lnlFno) ((0, J), 5) in the whole region of 5, where fidelity 
is independent of c, i.e. in the small-system regime, see Fig. [3] (bottom left). 

On the other hand, for given 5 and sufficiently large L the grid of wave numbers samples 
densely the vicinity of ir/2, where f(k) varies strongly and is sensitive to c. Consequently, 
— ln9 r (i i o)((0, J), 8)) is sensitive to c and scales with the system size as L. This is the 
macroscopic-system behavior. In Fig. |3l where L is fixed, this behavior is well visible for 
sufficiently large 8. 

The small-macroscopic system crossover is given by those 8 and L that satisfy the 
crossover condition ([2]), which on substituting ( 126|) assumes the form 



\8\L 



1 + J 2 



and for sufficiently large \J\ simplifies to 



\8\L 



1 J 



\J\ 



(34) 



c+1 



(35) 



Numerical calculations show that crossover condition ( 1351) is fairly well satisfied. 

In the macroscopic-system regime, we have been able to find analytical approximations 
to — ln3 r ( 10 )((0, J), 8)) for c > 1 and for c = 0. Specifically, for c > 1 



ln% )((0,J),5)) 



\_8\L_ 
Ac\J\ 



(36) 



and for c = 0, 



-lnJ (1 , 0) ((0,J),5)) 



-\5\L 



2tt(1 + J 2 ) 



(2 + J 2 )ln 1 



1 

J 2 



2\J\ I 7r — arctan 



J\ 



(37) 



For < c < 1 we have not been able to find approximate formulae analogous to the above 
ones. Formulae (136]) and ( J3TJ) compare well to numerical data shown in Fig. |3j Finally we 
note that whether it is a small- or macroscopic-system regime, fidelity depends only on the 
product 8L, see Fig. [31 bottom right panel. 

Thus, near A c =(0, J), with J away from 0, in transitions across the /x=0 critical line 
scaling laws ([6]) and ([7]), together with crossover condition ([2]), ( 135]) set by £ ( 126]) are obeyed. 

3.2 Across the /i=0 critical line: A c =(0, J), J close to 0, e=(l,0) 

In this section, we consider J being close to zero, say \J\ < 0.1, so the critical point 
A c = (0, J) is in a close vicinity of the multicritical point (0,0), and we consider only 
c > 0. As in the case considered in previous section 13.1] k c = n/2, f{k) is given by for- 
mula (130]) . but the character of f(k) changes significantly: the plots of f{k) acquire well-like 
shapes, see Fig. HI with the positions of the "wells" depending weakly on J. To reveal 
further properties of f(k) we make the corresponding plots of — In f(k) (right panel of Fig. 
H]). Away from tt/2 one can see, as in section [3~T1 — lnf(k) ~ k~ 2 behavior, which implies 
— ln3 r (i i o)((0, J), 5) ~ L 2 scaling, for sufficiently small 5. In the small-system regime, fidelity 
is described approximately by formula (133]) obtained in section |3~T1 the quality of this approx- 
imation is shown in Fig. [6l bottom left panel. We note that in distinction to the situation 
encountered in section I3TT1 the points of maximum curvature of the plots of — ln/(/c) (Fig. 
HI) depend weakly on J. Another view of the plot of — In f(k) for c > 1, with some cut-off of 
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Figure 3: (Color online). The plots of — ln3 r ( 10 )((0, 2), 8) versus 5 in doubly logarithmic 
scale, for two system sizes, L = 10 3 and L = 10 4 , and three values of c; for each system size, 
from bottom to top: c = 2 - magenta dotted line, c = 1 - green dashed line, c = - red line. 
In top right panel the data for L = 10 4 have been divided by 10. As a result, for sufficiently 
large 5 (i.e. in the macroscopic-system regime) the plot for L = 10 4 collapsed onto that for 
L = 10 3 . The blue dotted and black dotted straight lines represent analytical approximations 
(I37|) and ( )36|) . respectively, which agree well with numerical data in the macroscopic-system 
regime. In bottom left panel the data for L = 10 4 have been divided by 10 2 and for 
sufficiently small 5 (i.e. in the small-system regime) the plot for L = 10 collapsed onto that 
for L = 10 3 . Here, the blue dotted straight line represents analytical approximation ( 1331 . 
which again agrees well with numerical data in the small-system regime. In bottom right 
panel, in the plot for L = 10 4 variable 5 has been divided by 10 which resulted in collapse 
onto the plot for L = 10 3 . Scaling is consistent with v — 1, i.e. — ln3 r ( 10 )((0, J), 5)) ~ 5 2 L 2 
in small-system regime and — ln9 r ( lj0 )((0, J), 5)) ~ \5\L in macroscopic-system regime. 
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extremely small values, is shown in Fig. El for a sequence of parameter 5 values. For given L 
and sufficiently small S the width of the " dome" is much smaller than 2tt/L that separates 
consecutive wave numbers and its base does not incorporate any wave numbers (meaning 
that the contributions — In f(k) to — In 3 r ( 10 )((0, J), 5) are negligible). With increasing 5 the 
"dome" gets wider and "moves" towards the closest to it/2 wave number k\ = it/2 + tt/L. 
When the dome "passes" over kx, the contribution — ln/(7r/2 + tt/L) increases, reaches a 
maximum and decreases to a very small value. As a result, a first oscillation with pro- 
nounced maximum appears in the plot of — ln3 r ( 10 )((0, J), 5) for c = 5, see Fig. [61 Then, 
this scenario is repeated but for the next wave number k 2 = vr/2 + 3tt/L; consequently the 
second maximum in those plots develops. With increasing 6 the "dome" base gets wider 
than the separation of wave numbers, 2tt/L, so more than one wave number can fall under 
the " dome" . Therefore, the oscillations become less pronounced, their amplitude decreases, 
eventually a simple — ln9 r (i i0 )((0, J), 5) ~ 5 scaling results for sufficiently large S. 

We can easily determine the locations in 5 and the values of the described maxima of 
— In3 r ( 1]0 )((0, J), 5). For sufficiently small \5\, maxima of — ln/(/c) (in continuous variable 
k) are attained at the points 



n'i 



tt/2± 



c 2 -l 



J 2 



6, 



and their value is 



-In /(** 



In 



1 f 1 yV -!)(! + J 2 ) -c 



(3* 



(39) 



Therefore, for finite L the maxima of — ln3 r (io)((0, J), 5) are attained at the points S n 



± 



1 + J 2 (2n + 1)tt 



c 2 -l 



L 



n 



0,1, 



(40) 



see Fig. |6l bottom right panel. In the macroscopic-system regime, for J close to 0, we 
obtained analytical approximations to — ln3 r (io)((0, J), 5) for different c. For c = formula 
(l37l) of the previous section holds. Then, for c = 1 



ln% i0) ((0,J),5) 



L|*| 



7T 



1 4 5 t2 



and for c > 1 



-ln% )((0,J),5) 



L|tf| 



7T 



In 1 



J 2 c 2 



I Jlcarctan' 



Ulc 



-2 



(41) 



(42) 



The plots obtained from formulae ( 1371) . ( 14TT) and ( 1421) compare well with corresponding plots 
obtained numerically, see Fig. El top right panel. 

The plots in Fig. [6] make clear that the location of the small- macroscopic-system 
crossover region can be identified with the location of the first abrupt increase in the plot 
of — ln3 r ( 10 )((0, J), 5) versus 5. We have verified that this abrupt increase is determined by 
the position of the point of maximum curvature in the plots of — In f(k). For c > 1, it is the 
location of the first maximum in the plot of — In 3 r (i ) o)((0, J), that can be identified with the 
crossover region. Therefore, the crossover condition for the transition from (SL) 2 scaling to 
\5\L scaling of — ln3 r ( li0 )((0, J), 5) reads 



\S\L 



J 2 



-7T. 



(43) 
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Figure 4: (Color online). The plots of / (left panel) and — In / (right panel) versus tt/2 — k 
for three values of c: c = - red solid line, c = 1 - green dashed line, c = 5 - magenta dotted 
line. Right panel is in doubly logarithmic scale; the blue dotted straight line of slope —2 
indicates the —In f(k) ~ n~ 2 behavior, red and blue vertical lines indicate the location of 
the wave number 7r/2 — tt/L that is the closest to ir/2 for L = 600 and L = 100, respectively. 



For sufficiently small \J\, when the considered critical point is in close proximity to the 
multicritical point (0,0), the crossover condition (1431 is almost independent of J. Formula 
(1431) is in contrast with the expected condition, given by formula (1341) . which for sufficiently 
small \J\ reads: \8\L ~ l/(c + 1)|J|. This result can be interpreted as the appearance of 
a new, J-independent for sufficiently small \J\, characteristic length in the system, which 
we denote £'. The appearance of £' can naturally be attributed to the influence of the 
multicritical point. 

Thus, near A c =(0, J), with J close to 0, in transitions across the /i=0 critical line scaling 
laws (jBD and (j7|) are obeyed but crossover condition (J2J), (1331) set by £ fT2B1) is not. 



3.3 Along the J-axis: A c =(0, J), arbitrary J, e=(0, 1) 

Along the J-axis, that is a critical point and its neighborhood are contained in J-axis, we 
encounter a particularly simple situation. Since [i = jl — 0, f{k) simplifies to a constant: 



f(k) = I J i + l + [J+(c-l)S\[J+(c+l)S\ 

2 I V / (l + [J+(c-l)5] 2 )(l + [J+(c+l)5] 2 ) 

Expanding — In f(k) in 5 at zero we obtain 



(44) 



\n? m ((0,J),6) 



2 V 1 + J 2 



L + 0(S 3 )L. 



(45) 



Taking into account only the lowest order in 5 term, we obtain an accurate approximation to 
— In 3 r ( 0) i)((0, J), 5). The anomalous scaling (I45p holds in a wide range of values of parameter 
5. As shown in Fig. [7J this range covers the whole critical region, including the small and 
macroscopic system regimes, there is no small-system - macroscopic-system crossover. The 
plots of — ln3 r ( 0j i)((0, J), 5) versus 5, obtained for different values of J or L, after suitable 
rescaling according to (1451) . merge into one curve. 

Thus, in transitions along J-axis, near A c =(0, J), there is no small-system - macroscopic 
system crossover and scaling laws (jHJ) and ([7j) are not obeyed. 
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Figure 5: (Color online). For J = 10 2 , c = 5, L = 10 3 , the plot of —In/ versus 
tt/2 — fc (abscissa in logarithmic scale) for a sequence of 5; from left to right, 510 4 = 
5, 6, 8, 15, 20, 25, 37, 45, 50, 70, 100, 150, 200, 1000. For each 5 the plot is of steep-dome-like 
shape whose base width is approximately 25. The black points on abscissa mark the posi- 
tion of wave numbers, which are separated by 2tt/L. The black points on a dome indicate 
how many wave numbers are included in its base and what is their contribution to fidelity. 



3.4 The multicritical point A c =(0, 0), various directions e 

At the multicritical point the gap E\. closes at k c = ir/2. In a line neighborhood in direction 
e = (cos </>, sin 0), where < <ft < it/2, of the multicritical point the function f(k) assumes 
the form 



/(*) = 2 < l 



(cos k — /x) (cos k — p) + J J cos 2 k 



[(cos k — ji) 2 + J 2 cos 2 k] (cos k — fl) 2 + J 2 cos 2 k 



(46) 



where we set 



fj, = (c — 1)5 cos 4>, fi = (c + 1)5 cos i 



J = ( c - 1)5 sin 0, J= (c+l)5sin0, (47) 



with c > 0. 

If 5 is not too large, say 5 < 10 -2 , the plots of f(k) for different c are of well-like shape 
(see Fig. [S]). The wells are narrow and located close to k c = tt/2; their widths, calculated at 
the level f(k) = 1/2, are proportional to 5cos0, as a simple calculation reveals. Moreover, 
looking at the plots in Fig. [8] it is clear that for small systems, if there are no wave numbers 
inside the wells, f(k) equals approximately 1. Therefore, to get an approximate analytic 
expression for — In 3 e {\ c , 5) in the small-system regime we expand f(k) in 5 at zero, keeping 
only the lowest order term, which is /c-independent, 



l n 9 r (cos0,sin0)((O,O),5) 



sin 



-5 2 L. 



(48) 



This formula compares very well with numerical data shown in Figs. [HI HOI and HU Clearly, 
small-system scaling law ()6]) is not obeyed. 

For macroscopic systems, when for given 5 the system size L is large enough, it is the 
bottom of the well that plays a dominant role. In a narrow well this bottom is quite flat, so 
we expand f(k) in k at the minimum of the well, again keeping only the lowest order term. 
Then, the result, which is sensitive to the value of c, is integrated between the two points that 
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Figure 6: (Color online). Plots of — ln3 r ( lj0 )((0, 10~ 2 ), 5) versus 5 in doubly logarithmic 
scale, for two system sizes, L = 10 3 and L = 10 4 , and three values of c; for each system 
size, from left to right: c = 5 - magenta dotted line, c = 1 - green dashed line, c = - 
red line. In top right panel the data for L = 10 4 have been divided by 10. As a result, 
for sufficiently large 5 (i.e. in the macroscopic-system regime) the plots for L = 10 4 and 
for L = 10 3 merged into one curve. The straight lines in the upper part of the figure 
represent analytical approximations: formula ( 157)) for c = - blue dotted line, formula (14"T)) 
for c = 1 - red continuous line, and formula (142 p for c = 5 - black dotted line. All the three 
analytical approximations agree well with numerical data for sufficiently large S, i.e. in the 
macroscopic-system regime. In bottom left panel the data for L = 10 4 have been divided 
by 10 2 . As a result, for sufficiently small 8 (i.e. in the small-system regime) the plots for 
L = 10 4 and L = 10 3 merged into one curve. The blue dotted straight line represents 
analytical approximation ( I33p . which agrees well with numerical data in the small-system 
regime. In bottom right panel, variable 8 has been divided by 10. As a result, the plots 
for L = 10 4 and L = 10 3 merged into one curve. The blue dotted line indicates that the 
maxima of the oscillations in the plot for c = 5 are of the same height, given by ( 139)) . and 
are attained at points 5 n , given by ( HUjl . The orange vertical lines mark the positions of 
the maxima on the 5-axis. These positions are: 510 4 =6.3, 19.3, 32.1 (not shown), and 37. 
Scaling is consistent with v — 1, i.e. — ln3 r ( 10 )((0, J), S) ~ 5 2 L 2 in small-system regime and 
— In3 r (i )0 )((0, J), S) ~ \5\L in macroscopic-system regime. 
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Figure 7: (Color online). — In 3^0,1) ((0, J), 5) versus 5 in doubly logarithmic scale, for 
different c; from bottom to top: c = 5 - blue line, c = 1 - green line, c = - red line. 
The black dotted straight line of slope 2, coinciding everywhere with the red line and for 
In 5 < —4 also with green and blue lines, represents the first term of formula (14"5]) . Results 
consistent with infinite £. 
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Figure 8: (Color online) The plots of / (left panel) and — In / (right panel) versus k = ir/2—k 
for (j) = 60° and three values of c: c = - red solid line, c = 1 - green dotted line and, c = 5 
- blue dotted line. Right panel is in doubly logarithmic scale. 

are solutions to the equation f{k) = 1/2. For sufficiently small 5 the approximate solutions 
of this equation are given by the solutions k%,k2 of equations: cos/ci = jx and cosk 2 = ii. 
For all c, the result of integration is of the form C\8\ + C'\5\ In \5\, where C, C depend on c 
and 0. 

Specifically, for c < 1, the minimal value of f(k) is zero and is attained at n/2, and the 
lowest-order term of Taylor expansion at k = ir/2 is 



./' 



IX 



sin 



-AC 



.2 7 5 2 (l-c 2 ) 2 cos 4 , 
which after integration gives 

-L- 1 ln^ (cos0 , sin0) ((O,O),5)^C 1 |5|+C 2 |5|ln|,5|, 

with 

2cos</> / , 1 + c , . ,\ „ 2 cos i 



Ci 



TX 



2 — c In In sin < 

1-c 



Co 



TX 



(49) 



(50) 



(51) 
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Figure 9: (Color online) — In 3 r ( C os60 o ,sin60°)((0, 0), e ) versus L, in doubly logarithmic scale, 
for three values of c; from left to right: c = 5 - blue line, c = 1 - green line, and c = - 
red line. Two black dotted straight lines of slope 1 indicate that in the small system and 



macroscopic system regimes — In £F ( 



(cos 60°, sin 60°) 



60°) ((0, 0), e ) scales with the system size as L. 
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Figure 10: (Color online) Plots of — In 3 r ( C os60 o ,sin60°)((0, 0), 5) versus S in doubly logarithmic 
scale, for system size L = 10 6 and three values of c; from left to right: c = 5 - blue line, 
c = 1 - green line, and c = - red line. In left panel, the lower black dotted straight 
line of slope 2 is the plot of ( l4"8l) . and well agrees with numerical data in the small-system 
regime. The upper black dotted straight line of slope 1 indicates that in the macroscopic- 
system regime — ln3 r ( cos60 o iS i n 6o°)((0, 0), 5) does not scale as a power of S. Right panel shows 
the magnified view of the macroscopic-system regime. Here, the lower black dotted line is 
the analytical approximation (156]) to the blue line and the upper black dotted line is the 
analytical approximation ( !50l) to the red line. The uppermost magenta straight line of slope 
1 is a reference line. 
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Figure 11: (Color online) The system size is L = 10 6 . — In3 r ( cos8 9.9o,sin89.9°)((0, 0), 8) versus 
8, in doubly logarithmic scale, for three values of c; from left to right: c = 5 - blue line, 
c = 1 - green line, and c = - red line. The lower black dotted straight line of slope 2 is 
the plot of 



Then, for c = 1 



'(£- 



tan 



therefore 
with 



L-Mn?, 



(cos 0,sin^ 



l ((0,0),8)^C 3 \8\ + C 2 \8\ln\8\, 



2 cos <p 
C3 = ( 1 - - In sin < 



7T 



(52) 
(53) 

(54) 



Finally, for c > 1 we approximate the bottom of the well by a constant equal to the 
minimal value of the well, which for sufficiently small 8 is approximately equal to 



2x2-2 



min/(fc) « (cr-lprsin 



Then, integration leads to 

-L^lnS", 



with 



C, 



(cos 0, sin 0) 



2 cos 1 



71 



((0,0),8)^C 4 \8\+C 2 \8\ln\8\ 



In ((c 2 — 1) sin 



(55) 



(56) 



(57) 



In Fig. [TU] (right panel) it is shown that, for some range of 8 in the macroscopic-system 
regime, formula (|50|) well agrees with — ln3 r ( COS 6o°,sin60°)((0, 0), 8) numerically calculated for 
c = and c = 5. Similar results have been obtained for other values of angles </> (between 
and vr/2) and parameters c. However, if approaches tt/2, then in the macroscopic-system 



regime 



In J, 



(cos</>,sin< 



)((0, 0), 5) tends asymptotically, with increasing 8, towards the straight 



line describing small-system fidelity, given by (148 p . For very close to ir/2, this small-system 
fidelity is identical to that found at the line of critical points given by fi = 0, see Fig. [TTJ 

An inspection of figures Fig. [91 [TO], [TTJ shows clearly that the small- macroscopic-system 
crossover can readily be defined by the well-localized first "jump" in the plot of fidelity versus 
5 or versus L. Numerically we found that the position 8 of the "jump" scales with system 
size like L~ l , for a number of values of c. A simple argument provides an appropriate formula 
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- the crossover condition. For sufficiently small 5, there are no points of wave number grid 
inside the well in the plot of f(k), so — ln^cos^sm^CCO, 0), S) ~ 5 2 . When the first point of 
the wave number grid enters the well, — In 9 r (cos^,sin^>)((0, 0), S) jumps to a much greater value. 
This occurs approximately for those S and L that satisfy the equation cos(7r/2 — ac) = jl, 
with k = ir /L and jl = (c + 1)5 cos 0. That is, the crossover condition reads 



\S\L 



c+1 



7T- 



cos< 



(5? 



Formula ([58]) agrees well with our numerical data, in particular with those in Fig [T2l an d is 
consistent with v — 1. 

Summarizing scaling properties near multicritical point A c =(0, 0), we have found that 
neither small-system (jBJ) nor macrosopic-system (J7|) scaling laws are obeyed; the dependence 
of - In ^(cos^su^^O, ()),<!>) on \5\ is not power-like. 
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Figure 12: The location in 5 of the crossover region is linear in L~ l . The plot is made 
for c = and = 60°. Red crosses stand for numerical data, green dashed line represents 
formula ( 158|) . 



3.5 Across the J=0 critical line: A c =(/i, 0), 0<^-<l, e=(0, 1) 

In contrast to the critical points considered in previous sections, in this case the excitation 
gap E k closes at k c which depends on the system parameter fi; k c is the solution of the 
equation cos k c = fi. This fact has profound consequences for the behavior of fidelity. In a 
line neighborhood of a critical point (/i, 0), in direction (0, 1), the function f(k) assumes the 
form 



/(*) = 2 < l + 



(cos k — /i) 2 + J J cos 2 k 



[(cos k — fi) 2 + J 2 cos 2 k] (cos k — fi) 2 + J 2 cos 2 k 



(59) 



where we set 



J = (c-1)6, J=(c+l)S. 



(60) 



Plot of f{k) in a neighborhood of k c is shown in Fig. [131 left panel. For all c, the function 
f(k) is continuous and varies rapidly only in a neighborhood of k c , with the minimal value 
attaining zero if \c\ < 1; therefore in the small-system regime fidelity is insensitive to c. The 
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Figure 13: (Color online). The plots of / (left panel) and —In/ (right panel) versus 
(arccos/x — k) for three values of c: c = - red solid line, c = 1 - green dashed line, c = 2 
- blue dotted line. Right panel is in doubly logarithmic scale; the black-dotted straight line 
indicates the —In /(A;) ~ k~ 2 behavior. 
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Figure 14: (Color online). Left panel: plot of — In 3^0,1) ((1/3, 0), e -10 ) versus L in doubly 
logarithmic scale, for c = 1/2. Significant oscillations whose amplitude decreases with L 
for sufficiently large L are well visible. The black-dotted straight line of slope 2 represents 
formula (l6Tj) with adjusted numerical coefficient. The intersection of this line with the other 
black-dotted straight line of slope 1 indicates the crossover region between small-system and 
macroscopic-system behavior, which is manifested by the change of finite size scaling from 
L 2 to L. Right panel: plot of — m3v ,i) ((/•*, 0), e -10 ) versus L in doubly logarithmic scale 
for c = 1/2 and three values of \i\ from top to bottom: /i = 10 _1 - red line, \i = 10~ 2 - 
green line, fj, = 10~ 3 - blue line. There are two black dotted straight lines. One of them, 
whose slope is 1, given by formula (145]) with J = 0, coincides with the minima of the blue 
line for all L in the figure, and with minima of green and red lines for sufficiently small L. 
The other line, whose slope is 2, given by formula (16TT) with adjusted numerical coefficient, 
coincides with minima of the red line for sufficiently large L. The intersection of these two 
lines indicates, on increasing L, the crossover region between ultra-small-system regime with 
L-scaling and small-system regime with L 2 -scaling. 



20 



plot of — In f(k) in doubly logarithmic scale, Fig. [T31 right panel, reveals the — In f(k) ~ k~ 2 
behavior (k = arccos/i — k), sufficiently far away from arccos/i. However, this behavior does 
not imply — lnjF^o)^, 1), <5) ~ L 2 scaling in the small-system regime, as was the case in 
previous sections. 

The point is that for a typical critical point (/x, 0), when L is increased the grid of wave 
numbers "moves" with respect to k c . That is, if the points of the grid are initially, for some 
L, distributed symmetrically about k c , with increasing L their distribution becomes more 
and more asymmetric, and there is a wave number approaching k c arbitrarily close even for 
finite L. Further increase of L restores a symmetric distribution, and so on. If \c\ < 1, this 
"motion" of the wave- number grid with respect to k c causes strong variations of the values 
of factors f(k) in fidelity, which results in pronounced oscillations of fidelity, see Fig. HU 
In an analogous situation encountered in a one-dimensional XY model, similar oscillations 
have been ingeniously described by Rams and Damski [11] in a quantitative manner; they 
provided accurate approximate formulae for fidelity of sufficiently large systems. Therefore 
we do not dwell upon this point. 

For \c\ > 1, the oscillations of fidelity are much less pronounced, however for small \S\ 
and not too large systems they can still be significant. Fig. [T51 makes clear that for given S 
the oscillations are damped with increasing system size L, which can be seen in Fig. [14] 

As long as fidelity oscillates strongly with L, such notions as fidelity susceptibility and 
finite-size scaling have no meaning usually attributed to them. However, we can assign 
a meaning to those notions as follows. Fidelity attains a local maximum if k c is located 
approximately in the middle between two consecutive points of the grid, that is the clos- 
est to k c wave numbers are approximately k c ± n/L. These wave numbers give dominant 
contributions to fidelity for sufficiently small \S\. Taylor expanding hif(k c ± tt/L) in 5 and 
summing contributions from those two wave numbers we obtain the following approximation 
to fidelity: 

-ln^ m ((f,,0),5) ^ ^5Lf. (61) 

The above formula provides an envelope of minima of — ln9 r ( 0j i)((/x, 0), 5) as function of 
L. For given 5 and a range of system sizes, that are not too large and not too small (see 
our comment below), it fits very well numerical data, provided /i is sufficiently small, say 
ft < 0.1; for larger /i, the numerical coefficient 2/n 2 has to be increased suitably, see Fig. 
dU The set of S and L, for which formula (16"T|) approximates well the envelope of minima of 

— In 9 r ( 01 )((/i, 0), S) as function of L, for fixed S, can be qualified as the small-system regime. 
In Fig. [15] we show that, in the small-system regime, formula ([ST]) with adjusted numerical 
coefficient provides as well a good approximation to — ln3 r ( 01 )((/i, 0) as function of 5, for 
fixed L. Therefore, we can say, in the specified above sense, that in the small-system regime 

— In 3^0,1) ((/i, 0), 5) scales with S and L as (SL) 2 . 

In Fig. [TH (left panel) it is well seen that for fixed 5 the range of L, for which we deal with 
a small system in the above sense, is bounded from below. For smaller L, — In 3 r ( ,i)((/x, 0), S) 
scales as L, the system enters a new regime, which we call the ultra small-system regime. 
The crossover from the L 2 scaling to L scaling on decreasing sufficiently L, i.e. small- 
system-ultra small-system crossover, can be attributed to the proximity of the multicritical 
point. This is demonstrated in Fig. [HI (right panel) when the critical point (//, 0) ap- 
proaches the multicritical point (0, 0) the region of L scaling expands, incorporating larger 
values of L. In the limit /i — > the small-system regime (where L 2 scaling holds) disap- 
pears, lim^o — ln3 r ( 0jl )((/i, 0), 5) is well approximated by formula (J3HJ) with J = 0. Fig. [TJ] 
(right panel) indicates also that the small-system-ultra small-system crossover, on decreas- 
ing sufficiently L, occurs in a vicinity of the intersection of the two black dotted lines, given 
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Figure 15: (Color online). The system size is L = 10 4 . Plot of — In 3 r ( ,i)((2/3, 0), 8) versus 5 
in doubly logarithmic scale, for three values of c; from bottom to top: c = 5 - blue line, c = 1 
- green line, c = - red line. The three parallel black dotted straight lines of slope 1 indicate 
the macroscopic-system region, where — In 3^0,1) ((2/3, 0), 5) ~ 5. The black dotted straight 
line of slope 2 represents formula ( 161 j) with adjusted numerical coefficient and indicates the 
small-system regime, where it coincides with the three plots. 



analytically by formulae ( I45p with J = and (161 j) ; hence the considered crossover condition 
reads: 



fj?L 



(62) 



Concerning the macroscopic-system regime, — In jF( 0) i)((//, 0), 5) scales with 5 and L as 
\8\L, see Fig. [T5J Finally, we note that the small- macroscopic-system crossover takes place 
if the expected crossover condition L/£ e (A, S) ~ 1 holds. 

Summing up, in the generalized sense, explained above, in transitions across J=0 critical 
line, near A c =(/i,0), scaling laws (EJ) and (J7J) together with crossover condition (J2J) set by £ 
(I2"6"j) are obeyed with v = \. On decreasing sufficiently L, a crossover to a new regime, ultra 
small-system regime, takes place. 

3.6 Across the J=0 critical line, at the critical end point: A c =(l, 0), 
e=(0,l) 

As mentioned in section [2J in this case the critical index v — 1/2. At the critical point 
(1,0) the excitation gap E^ closes at k c — 0. In a line neighborhood of this critical point, in 
direction (0, 1), the function f(k) assumes the form 



/(*) 



where we set 



(cosfc — l) 2 + J J cos 2 k 



[(cos k — l) 2 + J 2 cos 2 k] 



J=(c-l)S, J 



(63) 



cos 



A;-l) 2 + J 2 cos 2 A; 



1)8. 



(64) 



Apparently, f(k) is continuous for all c, and away from k — it is insensitive to parameter 
c, see Fig. [T6j Moreover, away from k — it behaves like k~ A (right panel of Fig. [1~6|) . 
Consequently, in the small-system regime — In 3^0,1) ((1, 0), S) does not depend on c and 
scales with the system size as L 4 . To get an explicit approximation to — ln3 r ( 01 )((l, 0), S) 
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Figure 16: The plots of / (left panel) and — In / (right panel) versus k for three values of 
c; from bottom to top: c = - red continuous line, c = 1 - green dashed line, c = 2 - blue 
dotted line. Right panel is in doubly logarithmic scale and the black dotted line of slope —4 
indicates the region of —In f(k) ~ k~ 4 behavior. 

in the small-system regime, in the lowest order in 5, it is enough to use the approximation 
cosk ps 1 — k 2 /2, Taylor expand f(k) in 5 at S = 0, and take into account the contributions 
of the two closest to k c = wave numbers ±n/L: 



4 



ln%i)((l,0),«y)«-,^L 4 . 



7T J 



(65) 



In the macroscopic-system regime, the integral (128]) can be approximated by \S\ ' A(c), that 



is 



hx2 (M ((l,0),6)K i \6\ 1 ' 2 LA{ 



where the universal function A(c) is given by 



x 



A ^ vrio dX x* + 8cW 16(c 2 



(c 2 + V2c 2 - 1) 



i/i 



* - V2?3i)i/4 



8V2c 2 - 1 



(66) 



(67) 



Formula (|67|) holds for all c. The function A(c) is positive everywhere, at c = it amounts 
to cos(37r/8)/4, and at c = 1 its value is 2 _11//4 . We note that for c = 1, that is when one of 
the ground states used for calculating fidelity is the ground state at the critical point, A(c) 
is nondifferentiable. 

In Fig. [T7] we show that numerical data compare well with the approximations given by 
formulae ( 165]) and (166]) . Thus, in transitions across the J=0 critical line, near the critical 
end point scaling laws (Q and (JTJ) together with crossover condition (J2J) set by £ ( 126]) are 
obeyed with v = 1/2. 

4 Summary 

In this paper we have presented results of detailed studies of critical behavior of quantum 
fidelity in a BCS-like model. We have limited our presentation to a one- dimensional version 
of this model, where we have been able to derive a number of analytical results for cor- 
relation functions and quantum fidelity in critical regions. Needless to say that analytical 
results facilitate a lot understanding of numerical data. The presented results, obtained for 
a model that is out of the realm of paradigmatic quantum spin chains, can be thought of 
as complementary to those obtained by Rams and Damski [TU], [TT]. Results pertaining to 
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Figure 17: The system size is L = 10 4 . The plot of — ln3 r ( ,i)((l, 0) c , 5)) versus 5 for three 
values of c; from bottom to top: c = 2 - blue dotted line, c = 1 - green dashed line, c = - red 
continuous line. The black dotted straight line of slope 2 represents the small-system analytic 
approximation (|65|) . The two parallel black dotted straight lines of slope 1/2 represent 
the macroscopic-system analytic approximation f !66p . Scaling is consistent with v = 1/2, 



i.e. — ln3 r ( 01 )((l, 0), 5) ~ 5 2 L A in small-system regime and — ln3 r ( 0jl )((l, 0), 5) ~ l^l 1 / 2 // in 
macroscopic-system regime. 



Critical region 


Small s. 


Macroscopic s. 


Crossover 


V 


A c = (0,J), e = (1,0), large \J\ 


(SLY 


\S\L 


l ~ e i} 


1 


A c = (0, J), e = (1,0), small \J\ 


[5Lf 


\S\L 


l ~ e /(2) 


1 


A c = (0,J), e = (0,l) 


5 2 L 


5 2 L 


no crossover 


undef. 


A c = (0, 0), various e^ (0, 1) 


5 2 L 


\5\L + \5\\n\5\L 


L ~ ir/5cos<j) 


undef. 


A c = (//, 0), e = (0, 1), large |/i| (3) 


[SLf 


\S\L 


L~e 


1 


A C = (1,0), e = (0,l) 


5 2 L* 


\6\V 2 L 


L~e 


1/2 



Table 1: Resume of — lnjF^Ac, 5) scaling in small-system and macroscopic-system regimes, 
for critical regions investigated in the paper, except the case of A c = (/i, 0) being near the 
multicritical point (see section [375]) . W The correlation length £ is specified in sections |2] and 
13.11 ^ 2 > The length £' is defined in section 13.21 < - 3 - ) This scaling holds in the generalized sense 
(see section l3~5l) . 



the two-dimensional case, where almost no aid of analytical results is available, but still the 
exact diagonalizability of the model plays a key role, are in preparation and will be published 
elsewhere. 

We have verified critical scaling laws of quantum fidelity in the whole range of critical 
behavior: the critical scaling law (JSJ) in the small-system regime, where the effective lin- 
ear size of the system satisfies the condition L/£ e (\,5) ^C 1, then the crossover condition 
L/£ e (\,5) ~ 1, and finally the critical scaling law ((7j) in the macroscopic-system regime, 
where L/£ e (\,6) ^> 1. In the latter case, performing the thermodynamic limit, L — > oo 
for fixed 5, we have found that, for all the considered critical points, fidelity vanishes as 
— lnS^Ac,^) ~ L, that is the involved ground states become orthogonal; this is known as 
the Anderson orthogonality catastrophe [H] and has interesting consequences for various 
condensed-matter systems. 

First, in transitions across the J-axis [y = 1 - exact result), sections 13.11 and 13.21 and 
across the /i-axis but for the critical end point (A c — (1,0), v — 1/2 - numerical result), sec- 
tion l3.6[ fidelity scales according to the laws (jSJ) and ([7]), in the small-system and macroscopic- 
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system regimes, respectively. The behavior of fidelity in these transitions is fairly analogous 
to that in transitions across the line of critical points with specific value of transverse mag- 
netic field (g = 1) in anisotropic XY chain [TTJ. The crossover condition ([2j) with the 
correlation length £ given by (126]) is satisfied provided that the underlying critical point is 
not in a close proximity with the multicritical point. In the opposite case a characteristic 
length £', different from the correlation length £, determines the crossover condition. 

Second, we encounter a peculiar critical fidelity scaling in transitions that occur in 
neighborhoods of the multicritical point A c = (0,0), section [3T4l which is neither (J3D nor 
([7|). Specifically, in the small-system regime the anomalous scaling — ln3e(A c , 8) ~ 8 2 L 
holds. The crossover condition is consistent with v = 1. In the macroscopic-system regime, 
— hi.9r e (\.,8) is nonanalytic in 8 and the dependence on 8 is even not power-like. This 
behavior differs considerably from that in anisotropic XY chain [TT] . 

Third, an anomalous fidelity behavior takes place also in transitions along the J-axis, 
section 13.31 which is the line of critical points. Independently of the location, with respect 
to the multicritical point, of the critical point on this line, the anomalous scaling, found 
in the small-system regime in neigborhoods of the multicritical point, — In 3^ (Ac, 5) ~ 5 2 L, 
holds with no small- macroscopic-system crossover. The latter property is consistent with 
the correlation length being infinite. Such an anomalous scaling has been found also at the 
line g = 1 of critical points in anisotropic XY chain [11] , but only away from the multicritical 
point. 

Fourth, in transitions across the /^-axis (y — 1 - numerical result), section [331 fidelity 
as function of L may exhibit large-amplitude oscillations in the small-system regime, which 
makes such notions as fidelity susceptibility and finite-size scaling ill defined. One can assign 
a meaning to them in the sense of the envelope of minima of — In 3" e (A c , S) as function of 
L, and then, in the new sense, the scaling law (jSJ) is obeyed. The crossover condition (J2J) 
and macroscopic-system scaling (j7j) hold as well but in the new sense. Interestingly, an 
influence of the multicritical point can be seen when the considered critical point approaches 
the multicritical one (i.e. \x approaches zero), with L and 8 fixed and sufficiently small. 
Then, the small-system L 2 -scaling ([6]) changes to L-scaling, as it is in a neighborhood of 
the multicritical point, which we call the small-system - ultra small-system crossover; the 
crossover condition reads fi 2 L ~ 1. 

Finally, for reader's convenience we present in Table [T] a resume of our results concerning 
critical scaling of — ln3e(A c , 8) in small-system and macroscopic-system regimes, in various 
critical regions investigated in the paper. 
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